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Abstract 



. We study the equation of state of kaon-condensed matter including the effects of 

I temperature and trapped neutrinos. Several different field-theoretical models for the 

I nucleon-nucleon and kaon-nucleon interactions are considered. It is found that the order of 

■ the phase transition to a kaon-condensed phase, and whether or not Gibbs' rules for phase 

equilibrium can be satisfied in the case of a first order transition, depend sensitively on the 



. choice of the kaon-nucleon interaction. To avoid the anomalous high-density behavior of 



previous models for the kaon-nucleon interaction, a new functional form is developed. For 
all interactions considered, a first order phase transition is possible only for magnitudes 
of the kaon-nucleus optical potential > 100 MeV. The main effect of finite temperature, 
for any value of the lepton fraction, is to mute the effects of a first order transition, so 



P3 . that the thermodynamics becomes similar to that of a second order transition. Above a 

critical temperature, found to be at least 30-60 MeV depending upon the interaction, the 
first order transition disappears. The phase boundaries in baryon density versus lepton 



' number and baryon density versus temperature planes are delineated, which are useful in 

^ understanding the outcomes of protoneutron star simulations. We find that the thermal 

effects on the maximum gravitational mass of neutron stars are as important as the 
effects of trapped neutrinos, in contrast to previously studied cases in which the matter 
contained only nucleons or in which hyperons and/or quark matter were considered. Kaon- 
condensed equations of state permit the existence of metastable neutron stars, because 
the maximum mass of an initially hot, lepton-rich protoneutron star is greater than that 
of a cold, deleptonized neutron star. The large thermal effects imply that a metastable 
protoneutron star's collapse to a black hole could occur much later than in previously 
studied cases that allow metastable configurations. 

PACS numbers(s): 13.15.+g, 13.75.Jz, 26.60.+C, 97.60.Jd 
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I. INTRODUCTION 



It is believed that a neutron star begins its life as a proto-neutron star (PNS) in the aftermath of a 
supernova explosion. The evolution of the PNS depends upon the star's mass, composition, and equation 
of state (EOS), as well as the opacity of neutrinos in dense matter. Previous studies have shown 
that the PNS may become unstable as it emits neutrinos and deleptonizes, so that it collapses into a black 
hole. The instability occurs if the maximum mass that the equation of state (EOS) of lepton-rich, hot 
matter can support is greater than that of cold, deleptonized matter, and if the PNS mass lies in between 
these two values. The condition for metastability is satisfied if "exotic" matter, manifested in the form of 
a Bose condensate (of negatively charged pions or kaons) or negatively charged particles with strangeness 
content (hyperons or quarks), appears during the evolution of the PNS. 

Even if collapse to a black hole does not occur, the appearance of exotic matter might lead to a 
distinguishable feature in the PNS's neutrino signature {i.e., its neutrino light curve and neutrino energy 
spectrum) that is observable from current and planned terrestrial detectors. This was investigated recently 
by Pons et al. who studied the evolution of a PNS in the case where hyperons appeared in the star during 
the latter stages of deleptonization. Although the possibility of black hole formation was first discovered in 
the context of kaon condensation in neutron star matter [Q, a full dynamical calculation of a PNS evolution 
with consistent EOS and neutrino opacities in kaon condensed matter has not been performed so far. One 
of the objectives of this paper is to investigate K~ condensation in finite temperature matter, including 
the situation of trapped neutrinos in more detail. An impetus for this study is the recent suggestion that 
a mixed phase of kaon-condensed and normal matter might exist which could greatly affect the structure 
and its neutrino opacity 0. Another objective of our study is to identify differences in thermodynamic 
quantities such as the pressure, entropy or specific heat that might produce discriminating features in 
the star's neutrino emission. In separate works, we will examine neutrino interactions in kaon-condensed 
matter and neutrino signals from PNS evolution calculations in a consistent fashion. 

Since we wish to isolate the aforementioned effects due to kaons in this paper, we deliberately exclude 
consideration of hyperons. Hyperons and kaons were considered together in Refs. and [§]. Hyperons 
tend to delay the appearence of kaons in matter, especially if the S~ appears first. However, the E~ 
couplings are not as well determined as those of the A and even in this case the data are restricted 
to nuclear or subnuclear densities. Relatively small variations in the coupling constants can lead to a 
situation where the threshold density for the appearance of S~ particles is larger than that for kaons. 
These uncertainties remain unresolved; further hyper-nuclear experiments are needed to pin down their 
couplings. 

The original investigations of kaon condensation in neutron star matter {e.g. Refs. [^-11] and its 
astrophysical conseqences employed a chiral SU{?>)l x SU{?>)r model in which the kaon-nucleon 

interaction occurs directly via four point vertices. However, one can also employ an indirect, finite-range 
interaction which arises from the exchange of mesons. Several studies have been performed along these lines 
|^J^J^,|l^,0] . Ref. found that the chiral and meson exchange approaches give similar results provided 
that the kaon-nucleon couplings are chosen to yield similar optical potentials in nuclear matter. Allowing 
kaons to interact via the exchange of mesons has the advantage that it is more consistent with the Walecka- 
type effective field-theoretical models usually used to describe nuclear matter [|r^ . In most studies of kaon 
condensation it has been found that the transition to a phase in which kaons condense is second order for 
modest values of the kaon optical potential, Uk-, of order -100 MeV. For magnitudes of Uk well in excess 
of 100 MeV, however, the phase transition becomes first order in character. Even when the transition 
is first order, it is not always possible to satisfy Gibbs' criteria for thermal, chemical and mechanical 
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equilibrium, so a Maxwell construction, which satisfies only thermal and mechanical equilibrium, was 
sometimes employed to construct the pressure-density relation. 

Recently, Glendenning and Schaff'ner-Bielich (GS) ^ modified the meson exchange Lagrangian in 
such a way that the Gibbs criteria for thermal, chemical and mechanical equilibrium in a first order phase 
transition was possible. The extended mixed phase of kaon-condensed and normal matter which results 
produces a qualitative difference for the structure of a neutron star, since the EOS is softened over a wider 
region than in the case in which there is no mixed phase. This has implications for the mass-radius relation 
and the maximum mass, among other properties of the star. 

In this paper, we investigate the phase transition involving kaon-condensed matter and its influence 
upon the equation of state. We flnd that the precise form assumed for the scalar interactions (particularly, 
their density dependence), both for baryon-baryon and kaon-baryon interactions, determines whether or 
not the transition is flrst or second order, and, in the case of a flrst order phase transition, establishes 
whether or not a Gibbs construction is possible. Since the form of the scalar interactions is not experimen- 
tally well constrained at present, we have explored several different models in this study of the effects of 
kaon condensation on the EOS and the structure of a PNS. For each model, we have performed a detailed 
study of the thermal properties which are summarized in terms of phase diagrams in the density-lepton 
content and density-temperature planes. 

In Sec. II we present the various Lagrangians and derive exressions for the thermodynamic properties 
of each. We also develop the theoretical formalism necessary to describe baryons and kaon condensed 
matter in both the pure and mixed phases. This is followed by a discussion of the determination of the 
various coupling constants. Section III contains a comparison of the results for the EOS and for the 
structure of neutron stars for typical values of entropy and lepton content in a proto-neutron star as it 
evolves. Our conclusions and outlook for evolution of a proto-neutron simulation are presented in Sec. IV. 
In Appendix A, the extent of the correspondence between a meson exchange formalism and a chiral model 
to describe kaon condensation in matter is examined. The role of higher order kaon self-interactions in 
determining the order of the phase transition to a kaon condensed state is studied in Appendix B. 



II. THEORY 
A. Nucleons and Leptons 



We begin with the well-known relativistic fleld theory model of Walecka |]T5| supplemented by nonlin- 



ear scalar self-interactions |TB|. Here nucleons (n,p) interact via the exchange of a-, uj-, and p-mesons. 
Explicitly, the Lagrangian is 

= E ^ (^T'^^m - 9u.Yuj^ - gpYh^ .t-M*)N+ '^d.ad^a - \mla^ - U{a) 

-lu^^u^-' + im>^^^ - iB^, ■ B^"^ + imjb^ ■ , (1) 

where N is the nucleon fleld, the p-meson fleld is denoted by b^ and the quantity t is the isospin operator 
which acts on the nucleons. Scalar self-couplings ||T^, which improve the descripton of nuclear matter at 
the equilibrium density, are included in the potential U{a) = {bM / 3) {gaCrY + {c/ 4) {g^jcry , with M denoting 
the vacuum nucleon mass. The fleld strength tensors for the vector mesons are given by the expressions 
'^Aii/ = dfj,uj^ — dyUJ^ and 'B^y = d^h^ — d^h^. In the standard Walecka model the nucleon effective mass 

M*cj^ = M-g,a (2) 
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(the label GM refers to the Glendenning-Moszkowski parameters that we will use with this expression) . 
We shall also study an alternative form due to Zimanyi and Moszkowski (labelled by ZM) ||18|| : 



= E {(l + ^) ^ (^t'^^/. - g^r^^. - Qpl'^K NMn\ + \d^ad^a - Im'y - U{a) 

n,p ^ ^ ' 



-^u.^LU^'' + Imliu^LU^ - iB^, ■ B'^- + imjb^ ■ b'^ . (3) 

By redefining the nucleon field, — > {1 + g^^a/M)'^ N, the Lagrangian can be written exactly in the 
form Eq. (|^), but the nucleon effective mass becomes 

M*zM = M{l + g,a/M)-' . (4) 

For small values of a this is equivalent to the Walecka form. However the ZM effective mass has the 
property that, in the limit of large a, M^m remains positive whereas Mq^ can become negative ||7|,P^, 
which is unphysical. 

In the mean field approximation the thermodynamic potential per unit volume for both Lagrangians 

is 

^ = + U{a) - im^^^ - - 2T E / ^ In (l + e-^(^*-^^)) . (5) 



Here the inverse temperature is denoted hy (3 = 1/T, E* = yP + M*^ and the subscripts GM or ZM 
have been suppressed. The chemical potentials are given by 

IJ'P = J^p + Qlc^g + \gpbQ ; ^in = Vn + Qlc^g - \gpho . (6) 

Note that in a rotationally invariant system only the time components of the vector fields contribute to 
Eq. (^ and for the isovector field only the component contributes. The contribution of antinucleons is 
not significant for the thermodynamics of interest for a PNS and is ignored. 

Using VLj^, the thermodynamic quantities can be obtained in the standard way. The nucleon pressure 
is Pj\f = —Qj^/V, and the number density n„^p and the energy density ej\f are given by 



i=n,p 



(2vr)= 



where the Fermi distribution function frix) = (e^^ + 1) ^. The entropy density is then given by SV = 

The contribution from the leptons and antileptons is adequately given by its non-interacting form, 
since their interactions give negligible contributions Thus the thermodynamic potential per unit 

volume of the leptons and antileptons is: 

^ = - E r^?^ / [in (l + e-^^^-'^^)) + In (l + e-^(^^+'^^))] , (8) 



where and p^ denote the degeneracy and the chemical potential, respectively, of leptons of species £. 
The degeneracy gn is 2 for electrons and muons and it is 1 for neutrinos of a given species. Since the 



4 



star is in chemical equilibrium with respect to the weak processes p + i~ ^ n + ui, where the lepton i is 
either an electron or a muon, the chemical potentials obey fi^ — ji^^ = /ie — = IJ^n — IJ'p- If there are 
no neutrinos trapped in the star the neutrino chemical potentials fi^. are zero or, equivalently, the total 
neutrino concentration Y^, = 0, where we define the concentration for particle i to be Yi = ni/ijin + Up). 
The pressure, density and energy density of the leptons are obtained from Eq. in standard fashion. 



B. Kaons 



The two kaon Lagrangians of the meson-exchange type which have been previously suggested (in Refs. 
and PJ^, respectively), can both be written in the form 

Ck = d^K+d^K- - aK+K- + iX'^^K+d^R- - R-d^K^) , (9) 

where denote the charged kaon fields and we have defined the combined vector field 

Xf, = gujK^pi + gpKb^, ; (10) 

with ujfj, and 6^ denoting uj and the p° fields, respectively; g^iK and QpK are coupling constants. Since only 
the time components of the vector fields survive, in practice only Xq is non-zero. The two Lagrangians differ 
in the forms chosen for the quantity a. Both have the standard vacuum mass term, but the interaction 
terms differ. Specifically, Knorren, Prakash and Ellis (KPE) |^ take 

OlKPE = m*^;KPE = ^K- 9aKmKO- , (11) 

with rriK denoting the vacuum kaon mass, while Glendenning and Schaffner-Bielich (GS) choose 

«G5 = m*lGs - X^X^ = {rriK - y^KCxY - X^X^ . (12) 

Note that the coupling constant q^^k is defined here to be twice that defined in GS. A similar remark applies 
to the pNN coupling constant gp. It is remarkable, as pointed out in Appendix A, that to leading order 
in the kaon condensate intensity, the equations obtained with the chiral Kaplan-Nelson |^ Lagrangian at 
zero temperature agree precisely with those from the KPE Lagrangian. 

To see the significance of the term involving the vector fields in Eq. ([T^) , consider the invariance of 
the Lagrangian under the transformation K^e^^^. This allows the conserved kaon current density 

to be identified as 



= i{K+d^K- - K-d^K+) + 2X^K+K- . (13) 
Now for the combined GS Lagrangian, Cj^ + Ck, the equation of motion for the omega field is 

d'uj,^ + mlup = ^c. E - du^KJf. ■ (14) 



n.p 



Since the nucleon current Y.n,pX'~i^N is conserved, as is the kaon current J^, taking the divergence of Eq. 
(p!^) immediately yields d^uj^ = (and similarly for the p field). This is the required condition for a vector 
field so as to reduce the number of components from four to three. On the other hand, axpE does 
not contain an X^X^ term, so that the kaon current does not appear on the right hand side of Eq. (p 
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At the mean field level, however, where the vector fields are constants, any derivative is necessarily zero 
so that the divergence condition is automatically satisfied. 

For the coupling of the kaon fields to the scalar a field, KPE use a linear coupling, whereas GS have 
an additional quadratic term. There is little guidance on the form that should be used to generate the 
kaon effective mass so the choice is somewhat arbitrary, although, as we shall see, it can significantly affect 
the thermodynamics. Both the KPE and GS choices lead to problems for sufficiently large values of the a 
field; in one case the effective mass becomes imaginary, in the other it becomes negative. We are therefore 
led to consider a third form in the spirit of the ZM model for nucleons. For specificity we start with the GS 
Lagrangian which can be written = D ^K^ D^^* K~ — m*^.QgK^K~ in terms of a covariant derivative 

= + iXfj_, and replace it by 

C'k=(i + D^K+D^'*K- - mlK+K- . (15) 

We observe that the form of Cj^ above is one of many possibilities. Making the transformation 

(1 + ^g„K(^ I f^K)~^ and noting that a is a constant mean field, the kaon Lagrangian can be put in the 

form of Eq. (P) with 

aTW = m%Tw - ^M^^ = < (l + ' - X,X^ . (16) 

The label TW denotes "this work". While Eqs. (|l|), (0) and (|16D all give m*^ ^ rriK — \guKO' for small 
0", they differ at order and beyond, i.e., for large values of a. 

The kaon partition function at finite temperature can be obtained for a Lagrangian of the form (P) by 
generalizing the procedure outlined in Kapusta First, we transform to real fields (pi and 02, 



ir± = (01±202)/V^, (17) 

and determine the conjugate momenta 

TTl = <9o0l - Xo02 ; T^2 = <9o02 + ^001 • (18) 

The Hamiltonian density is Tix = tti^q^i + 7r25o02 — ^k, and the partition function of the grand canonical 
ensemble can then be written as the functional integral 

= J [d7Ti][d'K2] J ^ [d(f)l][d(j)2\ exp (^^^^ ~ ^^'^ ^ fJ'M'^i^ (t>i) 

(19) 

Here the fields obey periodic boundary conditions in the imaginary time r = it, namely 0i(x, 0) = 0j(x, /?), 
where (5 = 1/T. The chemical potential associated with the conserved kaon charge density is denoted by 
fx and chemical equilibrium in the reaction e" ^ K~ + requires that 

The Gaussian integral over momenta in Eq. (^) is easily performed. Next the fields are Fourier 
decomposed according to 

0i = /^cosC + J^Ee*^^"''""^Vi,n(p) ; 02 = /^sinC + J^Ee^^^"^""'V2,n(p), (20) 

' n.p ' ra,p 
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where the first term describes the condensate, so that in the second term 0i,o(p = 0) = 02,o(p = 0) = 0. 
The pion decay constant / has been inserted so that the condensate angle 9 is dimensionless. The 
Matsubara frequency uJn = 27rnT. The partition function can then be written 



Zk = N f l[[d(pi,n{p)][dhAp)](i'^ , where 

n,p 



n,p 



l-nl-pj,02, 



2{n + Xo)uj„ 



-p) p 



V(P) 



-2(/i + Xo)uJn oj^+p^ + a - 2/iXo - /x^ 

is a normahzation constant. We define the K"^ energies according to 



(21) 



(22) 



so that the three approaches give 



^KPEi.P) = + T^K;KPE + ^0 ± ^0 

'^Gs(p) = \Ip^ + ^K.GS ± ^0 



(23) 



Using the definition (^) and suppressing the explicit dependence of cu^ on p, the determinant of D is 



givmg 



V 



det D = ujI + {uj- - /i)2 < + (cu+ + /x)^ 



1 '7 1 

^ - |(/^)'(a-2/iXo-/i2) + — -^Indet D, 



2/?V^ 



(24) 



(25) 



n,p 



where the normalization constant N has been dropped since it is irrelevant to the thermodynamics. Per- 
forming the sum over n and neglecting the zero-point contribution, which contributes only beyond the 
mean field approach and in any case is small p2[, we obtain the grand potential for the kaon sector: 



-i(/^)2(a-2/iXo-/i2)+T 



(2^ 



ln(l -e 



+ ln(l - e 



-/3(a;++At) 



(26) 



The kaon pressure, Pk = —^k/V, and the kaon number density is easily found to be 
UK = {fefifx + Xo) + where = / ^slfBi^' " " /^^(^^ + 

and the Bose occupation probability /b(x) = {e^^ — 1)^^. The kaon energy density is 
SK = \Ue)\a + ^,[uJ~ip)fBiuJ- - /.) + i^^{p)fBico^ + /x)] , 

and the kaon entropy density is Sk = (3{sk + Pk ~ (^''^k)- 



(27) 



(2J 
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C. Equations of Motion 



It is useful first to define the quantity 

= / (p' + « + ^o)'" [/^(^" - + /^(^^ + • (29) 

Then the mean cj, p and a fields, as well as the condensate amplitude 9, determined by extremizing the 
total grand potential fitotai = + + can be written 



m>o = 9.{n, + n„) - g^j, { KfOY + n"/ - X.A"/ - \[{fef + ^F]^ 



= \9p{np - nn) - 9pK + - X,Al^ - \[{fdf + 



J 

So; 1 
da 



6{fi^ + 2fiXo -a) = e[fi - cu"(0)][/i + c<;+(0)] = . (30) 

The derivative da/dXo is zero for the KPE case and —2Xo for the GS and TW cases. The derivatives 
with respect to the a field are 

da ' da ^"K M 

daxpE _ dacs _ * darw _ i^K;Twf /oin 

da da ' da mj^ 

Note that the last of Eqs. (^) yields either 6 = (no condensate) or the condition for a condensate to 
exist. Since fi is positive here, we only have the possibility of a K~ condensate with fi = uj~{0). Note also 
that the contribution of the condensate to the kaon pressure Pk vanishes, as it should. 

The remaining condition to be imposed is that neutron star matter must be charge neutral. For a 
single phase this implies 

Hp-riK -ne-nf, = 0, (32) 

where Ue and are the net negative lepton number densities. The mixed phase thermodynamics is 
discussed below. 

The sum of the nucleon and kaon energy densities can be simplified somewhat by using the equations 
of motion. This gives 



0^ + Imlbl + 2 E / -0^,E*ME* - u^,,) + UOf {a - \X, 



da 

dX'a 



n,p 

+Xonl^ - XoAF (xo + + / ^^[uj-{p)fB{u;- -p)+ u^Hp)fB{u^^ + /.)] • (33) 
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D. Mixed Phase Thermodynamics 



In the theory discussed above there are two independent chemical potentials, which we can take to be /i„ 
and n, each connected with a conserved charge, the baryon number and charge of the system, respectively. 
Glendenning p3| pointed out that in the presence of a first order phase transition, conservation laws must 
be globally, not locally, imposed, if possible, in the mixed phase region. A Maxwell construction would have 
been appropriate had there been just a single conserved charge. However, relaxing the condition of local 
charge neutrality does not guarantee that the model Lagrangian, solved in the mean field approximation, 
will provide a description of the mixed phase, which is only possible if the Gibbs criteria can be satisifed. A 
simple, yet general, procedure to check if the Gibbs criteria can be fulfilled by a specific model is discussed 
in Appendix B. 

Denoting the phase containing a condensate with a subscript 6, and the phase without a condensate 
with 6 = 0, the total pressures in the two phases must be equal 

-P0=o(yUn, /i, T) = Peifin, Ai, T) . (34) 

Each of the chemical potentials is the same in the two phases. If the volume fraction of the non-condensed 
phase is then global conservation of charge requires 

X[np -riK-rie- r2^]e=o + (1 - x) hp - rix - - = . (35) 

The densities of the individual species in the mixed phase are evident here. The total energy density is 
the weighted sum of the two phases 

£ = X£e=oil^n, yu, T) + (1 - x)£^e(/in, /U, T) . (36) 
The total entropy density is obtained through a similar equation. 



E. CoupUng Constants 

In the effective Lagrangian approach adopted here, knowledge of two distinct sets of coupling constants, 
one parametrizing the nucleon-nucleon interactions, and one parametrizing the kaon-nucleon interactions, 
are required for numerical computations. These are associated with the exchange of a, uj and p mesons. 
We consider each of these in turn. 



1. Nucleon Couplings 

The nucleon-meson coupling constants are determined by adjusting them to reproduce the properties 
of equilibrium nuclear matter at T = 0. The properties used are the saturation density, no, the binding 
energy/particle, Ea, the symmetry energy coefficient, agymi the compression modulus, K, and the Dirac 
effective mass at saturation, M* . Not all of these quantities are precisely known and the values we choose 
are listed in Table I. For completeness, we list the equations needed to obtain the coupling constants, 
assuming that the scalar self-coupling has the form ?/($) = (6M/3)$^ + (c/4)$^, where $ = g^a. From 
the equation of motion for the ujq field and the fact that the pressure is zero at saturation density in nuclear 
matter, the value of g^jm^^ is given by 
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M -E 



A 

no 



EX 



(37) 



where E*p = ykp + M*^ and no = 2A;|,/(37r^). The p meson couphng constant can be determined for a 
given symmetry energy through the relation 



9l 
mi 



no 



2a 



sym 3^^^ 



(38) 



An expression involving the compression modulus can be deduced by differentiating the a equation of 
motion: 



'9^ 



9noM 



*2 



Ef[K + %Ea + E*p- M)] - 3klE*p 



no 
E*p 



M* 



df rf2f/($o) 



'd<^o 



(39) 



Here $0, the value of $ at saturation density, is obtained directly from the Dirac effective mass. The 
function /($o) = —dM*/d^o depends on the particular expression used for the effective mass. The scalar 
density n^ = {M* /7i'^){kFEp — M*'^\n[{kp + Ep)/M*]}. The a equation of motion at saturation can be 
written in the form 
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no 



E*p + EA-M)+ns (iM* + <l>o/ 



(40) 



which together with Eq. (|39D allows the a coupling to be obtained. Finally the constants appearing in the 
scalar self-coupling t/($) are determined from: 



1 



'2 



d$2 



2M$ 
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rf2f/($ 



d<l>l 



3c<l>l 



(41) 



The constants determined in this way are given in Table I. Note that in principle the potential should be 
bounded from below for large values of the a field requiring c to be positive; this is the case for the ZM 
model. 



2. Kaon Couplings 

In order to investigate the effect of a kaon condensate on the EOS in high-density matter, the kaon- 
meson coupling constants have to be specified. Empirically known quantities can be used to determine 
these constants, but it should be borne in mind that laboratory experiments give information only about 
kaon-nucleon interaction in free space or in nearly isospin symmetric nuclear matter. On the other hand, 
the physical setting in this work is matter in the dense interiors of neutron stars which has a different 
composition and spans a wide range of densities (up to ~ 8no). Therefore, kaon- meson couplings as 
determined from experiments might not be appropriate to describe the kaon-nucleon interaction in neutron 
star matter, and the particular choices of coupling constants should be regarded as parameters that have 
a range of uncertainty. 
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With the above caveats in mind, we now examine the relationship between the optical potential of 
a single kaon in infinite nuclear matter and the kaon-meson couplings in our Lagrangian. Lagrange's 
equation for an s-wave K~ with a time dependence K~ = k~{x.) e~*'^*, where E is the asymptotic energy, 
defines the optical potential for our Lagrangian (H) according to 



[V^ + E^- m\] k-{x) = [-2XoE + a - m\] k~{x) 

= 2 niK Uk k~ (x) . (42) 

In nuclear matter, feg = 0, so for a kaon with zero momentum {E = rriK) the optical potential is 

a — m\ , , , 

Uk = -i: - QuKOOo ■ 43 

Utilizing the functional forms for a in Eqs. ([TT|), (p!2D, and (p!6D, the optical potentials for the KPE, GS 
and TW models are easily obtained. For the KPE case this may be written exactly as 

Uk^^ = -\g.K(J - QukuJo , (44) 

whereas for the GS and TW cases there are higher order corrections in addition to the terms linear in the 
fields. We choose g^x to be gui/3 and gpK to be gp/2 on the basis of simple quark and isospin counting 
arguments. Note that this value for g^^x is also suggested by comparison to the chiral approach (see Ref. 
10] and Appendix A) and it leads to a —48.5 MeV contribution to the optical potential. The total optical 
potential is shown in Table II for various choices of the a coupling. The linear form Eq. (0), exact for 
KPE, is an accurate fit to the the GS and TW cases for moderate values of the optical potential. For 
orientation, chiral models suggest that the magnitude of the optical potential is at most 120 MeV [|^, 
while fits to kaonic atom data have been reported with values in the range 50-200 MeV |p5|-p^. We note 
that Glendenning and Schaffner-Bielich label their results according to values of the optical potential 
obtained in the linear approximation (henceforth, U^k^). In order to make an apposite comparison with 
their results, we will parametrize the kaon coupling for each model simply by specifying the value of 

rrlin _ ttKPE 



III. RESULTS AND DISCUSSION 

A. Zero temperature case 

The effects of kaon condensation on the EOS are more pronounced at zero temperature than at 
finite temperature, since the fraction of thermally excited kaons increases with temperature relative to 
the fraction of kaons residing in the condensate. We therefore begin by examining results for the zero 
temperature case. We have considered two different nucleon Lagrangians, GM and ZM, and three different 
kaon Lagrangians, KPE, GS and TW. Below densities of about 0.5no, matter is composed of neutron-rich 
nuclei immersed in a neutron sea. For this regime, we use the potential model results of Negele and 



Vautherin in the range 0.001 < n < 0.08 fm~ and those of Baym, Bethe, and Pethick for 
n < 0.001 fm~^. For cold stars, the EOS in this regime has little effect on maximum masses or stellar 
radii. Furthermore, since the entropy in the stellar mantle [n < no) is quickly radiated away in neutrinos, 
the EOS in this regime does not substantially affect the results of this paper. 

In Fig. 0, we compare the pressures for the different nucleon and kaon Lagrangians as a function of 
baryon density, Ub = rin + Up. The solid lines show results for both the pure nucleon and kaon condensed 
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phases with no attempt to enforce the Gibbs conditions of chemical and mechanical equilibrium. In all 
cases, a first order phase transition is found to occur, as long as the magnitude of the optical potential 
l^iPl — \9(tkO' + QuK^ is in excess of 100 MeV. Where possible, the pressure in the mixed phase obtained 
by imposing Gibbs' criteria for mechanical and chemical equilibrium is shown as a dashed line. For 
the GM+KPE, ZM+KPE and ZM+GS models it was not possible to satisfy Gibb's criteria, despite the 
occurrence of a first order phase transition for large enough jt/j^'^l. The reason for this is connected with 
the form of the kaon Lagrangian, as discussed below. We also point out in Appendix B that non-linear 
kaon self-interactions lead to a second order, rather than a first order, transition. 

The qualitative similarity of the results shown in Fig. |I] for the different nuclear Lagrangians enables us 
to simplify our analysis by allowing us to focus on three, rather than six, possible Lagrangian combinations. 
For a given kaon Lagrangian, fairly similar results can be obtained with different nuclear Lagrangians by 
making relatively small shifts in the kaon optical potential U^k^. The following discussion will therefore 
focus on the three cases GM+KPE, GM+GS and ZM+TW. The case GM+KPE is chosen to compare 
with the results of KPE, the case GM+GS is chosen to compare with the results of GS, and the case 
ZM+TW demonstrates the usefulness of Lagrangians in which anomalous values of effective masses are 
implicitly eliminated. The results for the model GM+GS shown here and elsewehere in this paper are 
identical to those found by GS for the same interactions. Note that in all models considered, the phase 
transition is second order in nature for moderately low values of the optical potential. 

In Figs. and |, the density dependence of the scalar, vector, and iso-vector fields, the electron 
chemical potential /Xg = yU = /i„ — /ip, the condensate amplitude 6*, and the nucleon and kaon effective 
masses are displayed in the pure nucleon and kaon condensed phases, ignoring any possible mixed phase 
for the present. For the optical potential chosen, = —120 MeV, a first order phase transition occurs 
in all three cases. After the onset of condensation a rapid change in the behavior of the electron chemical 
potential and some of the field strengths is seen to occur. The differences in the variation of the scalar 
{.QaK'^) and isovector {—gpRbo) fields between the models are particularly illuminating. For GM+KPE, the 
scalar field exhibits a relatively rapid increase with density after the onset of condensation. This in turn 
causes both the nucleon and kaon effective masses to drop rapidly with density. In fact, for sufficiently 
large density, the GM+KPE kaon effective mass vanishes (see Fig. ^). The variations of the effective 
masses in the models GM+GS and ZM+TW are more moderate. The variation of the isovector field 
with density, which in large part controls the variation of the electron chemical potential n and hence the 
electron concentration, is also more dramatic in the case of GM+KPE than in the GM+GS and ZM+TW 
models. Notice that in the KPE model it goes to zero for asymptotic densities (this follows from Eq. (|30D), 
so that the proton and neutron abundances become equal. This does not occur for the other two cases 
considered here. Finally, it is worth noting that in all three models the condensate amplitude rises rapidly 
once the threshold density is reached. 

We turn now to a discussion of the results obtained by imposing Gibbs' criteria for mechanical and 
chemical equilibrium at zero temperature. In Fig. ^, we show the chemical potentials associated with the 
two conserved charges, charge and baryon number, as functions of each other, for the model GM+GS for 
a kaon optical potential of t/j^ = —120 MeV. Quantities associated with the pure nucleon phase. Phase I, 
are shown as solid lines here and in subsequent figures. Phase II refers to the high-density phase in which 
nucleons and the kaon condensate are in equilibrium, and quantities associated with it are shown as dashed 
lines. Both phases, I and II, coexist in the mixed-phase region which is displayed as a dotted line. This 
figure illustrates the way a mixed phase is built from the two pure phases. For electron chemical potentials 
below the solid curve, matter is positively charged in phase I. A similar interpretation of positive or negative 
charge for below or above the dashed curve is not possible, since two different types of particles, kaons 
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and leptons, can furnish charge. In other words, a decrease in /i, or, equivalently, the number of electrons, 
does not necessarily lead to a positive net charge in phase II. For /i„ < 1165 MeV, only phase I with 
nucleons and leptons are present. For /i„ < 1310 MeV, a mixed phase of positively charged phase I and 
negatively charged phase II obeying the Gibbs' conditions (|34D is favored. Qualitatively, a similar situation 
is encountered in the construction of the mixed phase for the ZM+TW model, but the mixed phase region 
is quite small. As noted earlier, however, it was not possible to satisfy Gibbs' criteria for models with the 
kaon Lagrangian KPE. 

In Fig. 1^ we show the individual charge densities of phase I and II in the mixed phase, as a function 
of baryon density. The dotted curve in this figure shows the volume fraction of phase I. The results are for 
the GM+GS model with f/|? = -120 MeV (upper panel) and for the ZM+TW model with U^^ = -140 
MeV (lower panel). Near the lower threshold, matter in phase I is very slightly positively charged and 
occupies most of the volume. As the density increases, the volume fraction of phase I, X; decreases and 
its charge density increases. Note that the negative charge density of matter in phase II at the lower 
transition point, ~ 0.5 fm~'^, and the positive charge density of matter in phase I at the higher transition 
point, ~ 1 fm~^, are rather large in the case of GM+GS compared to the case ZM+TW. This is due to 
the stronger density dependence of the scalar and isovector densities in the former case. Note also that a 
first order transition allows for the existence of a very dense and nearly isospin symmetric matter in the 
mixed phase. 

In Figs. ^ and |^, we show the magnitudes of the various fields, the electron chemical potential, the 
nucleon and kaon effective masses, and the condensate amplitude for the GM+GS and ZM+TW models, 
respectively. Both models show the same qualitative behavior. At the lower phase boundary, in which 
phase II just begins to appear, the scalar field in phase II is much larger than in phase I and the condensate 
amplitude 6 in phase II takes a large value which decreases with increasing ub through the mixed-phase 
region. Thus, the effective masses of both kaons and nucleons in phase II are much smaller than in phase I. 
The densities demarking the mixed phase region and its overall extent are dependent upon the interaction 
models, and upon the assumed values of the kaon optical potentials, here taken to be —140 MeV in the 
case of GM+GS and —160 MeV in the case of ZM+TW. The region in density over which the mixed 
phase extends is much smaller in the latter case, chiefly due to the more moderate behavior of the scalar 
interaction with density variations in this case. 

It is instructive to compare the behavior of the two models at the threshold of the mixed phase region. 
Phase I will have a net small positive charge and a volume proportion x close to 1 (see Fig. This has to 
be counterbalanced by a large net negative charge in phase II since it is weighted by the small proportion 
(1 — x). Focusing on phase II, the condensate condition for the models GM+GS and ZM+TW from the 
last of Eqs. ( |5(]| ) is 



In order to ensure that uk > 0, the quantity (/i + Xq), and hence m*^, has to be positive definite. In 
the ZM+TW model the kaon effective mass is relatively large so that Xq is positive and therefore 6 is 
relatively small. On the other hand in the GM+GS model m}^ is quite small so that Xq is negative and 6 
has to be large. The negative value of Xq = Qujk^q + QpRbo implies a large negative value of 



/i + Xo = m;^ 

and the kaon number density, which has to be large, is 

nK = UOf{li + XQ). 



(45) 



(46) 




(47) 
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which is clearly sensitive to the value of gpK- In fact, if this coupling is reduced by more than about 
15% from our chosen value it is no longer possible to satisfy the Gibbs criteria. By comparing the pure 
phase results in Figs. ^ and ^ with the mixed phase results of Figs. ^ and 0, it is clear that substantial 
modifications of the various fields are required to satisfy Gibbs' criteria. 

We examine now the KPE model for which it was not possible to satifsfy the Gibbs' criteria. In this 
case, Eq. (|ll]) and the last of Eqs. (pOD leads to the condensate condition 

+ 2X0) = {m*j,f , (48) 

whereas the functional form of the number density of kaons is identical to that in Eq. (|^). Eq. (^8]) 
differs in important ways from Eq. (|45D . For the KPE model, even if /x + 2Xq is positive, has the 
proclivity to turn negative for large /i„ (or equivalently, for large baryon densities), leading to (m|^)^ < 
or imaginary values of the kaon effective mass m|^. This may be seen in Fig. ^ where we show the electron 
chemical potential as a function of the (negative) charge density in pure phase II for a typical value 
of the neutron chemical potential /i„ = 1250 MeV. It is now possible to understand qualitatively why a 
mixed phase cannot occur in the case of the kaon Lagrangian KPE. In comparison with the GM+GS and 
ZM+TW models, a distinctive feature of the KPE model is that fi decreases rapidly with the (negative) 
charge density. In constructing a mixed phase, we are attempting to balance the positive charge in phase 
I with the negative charge in the dense phase II in which the electron chemical potential, and hence the 
charge content in leptons, is rapidly decreasing towards zero. The balance never occurs, hence the failure 
to meet the Gibbs' criteria. In terms of compositions, the GS or TW Lagrangians introduce negative 
charges in matter by increasing the number density of kaons, while keeping the electron density nearly 
constant or even slightly increasing with the charge density. The KPE Lagrangian, however, rapidly 
substitutes electrons by kaons, which is detrimental to meeting the Gibbs' criteria. For these reasons, we 
will concentrate on results with the other two kaon Lagrangians in the remainder of this paper. 

The influence of the condensate on neutron star structure (at zero temperature) is shown in Fig. |^ 
in which the gravitational mass is displayed as a function of the star's central baryon number density 
(left panel) and its radius (right panel). For the models shown, the transition is first order and Gibbs 
equations for mechanical and chemical equilibrium are utilized. For all cases shown the central densities 
of the maximum mass stars lie in the mixed phase. The effects of the condensate are more evident in 
the case of the GM+GS model in which the mixed phase occurs over a wider region of density than in 
the ZM+TW model. When the effects of the softening induced by the occurrence of the condensate are 
large, the limiting mass and the radius at the limiting mass are reduced significantly from their values 
when the condensate is absent. Note, however, that the softening effects are limited by the constraint 
that the maximum mass must exceed that of the binary pulsar PSR 1913+16, 1.442 Mq. In the case of 
GM+GS, this constraint limits \U]^\ to be smaller than about 125 MeV. In such a case, the minimum 
radius achieved is not as small as in the case [/^" = —140 MeV, as shown in Fig. ^. The radii of stars 
with masses less than 1.2 Mq are not affected by the choice of the kaon Lagrangian or the kaon optical 
potential, since the condensation threshold is not reached in these cases. 



B. Comparison with other works 



The density dependence of m*j^/mK, —Uk and ujk have been investigated in other works [PB|- P5| , PT| , |52 
but for the most part either in isospin symmetric nuclear matter or pure neutron matter. In general, our 
results for mfj^/rriK with — f/]^'^ = 80 MeV are consistent with those of Refs. P^J23] (for an appropriate 
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comparison, our results are to be compared with results obtained without in-medium pion contributions 
in Ref. [^) and those of Ref. for nuclear matter at both n/riQ = 1 and 3. There is a relatively small 
change produced in going from nuclear matter to beta-equilibrated neutron star matter to pure neutron 
matter for the quantities m*p^/mK and —Uk- Note that a direct comparison of the real parts of the optical 
potentials between different calculations must also account for the fact that in obtaining fits to data, the 
imaginary parts are often found to be as large as the real parts, which indicates fragmentation of strength 
in the quasi-particle spectral function. 

Relatively larger variations are found in the kaon energies in matter with varying amounts of isospin as 
can be seen from Fig. ITUI. In this figure, the top panel provides a comparison of results for beta-equilibrated 



neutron-star matter for the GM+KPE, GM+GS, and ZM+TW models, respectively, for values of —U^^^ 
at the extreme ends considered here, namely, 80 and 120 MeV. The bottom panel shows results for the 
ZM-I-TW model for — f/^^ = 80 MeV, for pure neutron matter, neutron-star matter, and isospin symmetric 
nuclear matter, respectively. At nuclear density where the models are calibrated, u decreases by about 
a few MeV in going from pure neutron matter to neutron star matter and by about a few tens of MeV 
in going from neutron star matter to nuclear matter. With increasing density, these differences become 
progressively larger. This trend is chiefly due to the behavior of the vector fields in matter with different 
amounts of isospin. 

At this time, our results for the density dependence of uj can be compared with those of the potential 
models in Refs. [31 ,32| . For values of —U^k^ near the lower end of the range we explored, in the neighborhood 
of 80 MeV, the behavior of uj, for example, is quite similar to the potential model results. As the authors in 
Refs. [pl| , p2| indicated, kaon condensation may be unlikely in this case. However, the relevant comparision 
must also include the electron chemical potential /ig, since the density where u = He determines the 
onset of kaon condensation. As demonstrated in Ref. |3^, the behavior of /le for neutron star matter 
at high densities is determined by the density dependence of the nuclear symmetry energy (see also a 
similar discussion in Ref. |^). Potential model calculations (see, for example Ref. [Q) tend to have a 
relatively weak density dependence of the symmetry energy, which generally results in an onset of kaon 
condensation that is at a rather large density. In field-theoretical and Dirac-Brueckner-Hartree-Fock 153 



models, however, the symmetry energy varies relatively rapidly with density. These lead to smaller densities 
where kaon condensation occurs, for a given behavior of the kaon energy u. Furthermore, the calculations 
of Ref. 1^ have been performed only for pure neutron matter which further enhaces the values of u and 
discourages kaon condensation. In addition, as \U^j^\ is increased in magnitude in field-theoretical models, 
the role of kaons increases and u becomes progressively smaller as a function of density. Nevertheless, the 
lack of effective constraints at high density preclude choosing any model over another at this time. 

In summary, choosing values of —U^^^ near the lower end of the range we explored either lead to a 
second order phase transition or no transition at all in a neutron star, in which case the gross properties 
of the star are relatively unaffected from the case without kaons. On the other hand, values near the 
higher end of this range lead to a first order phase transition at a relatively low density, depending on the 
form of the interaction chosen, and a more pronounced effect on the star. Our aim has been to provide 
benchmark calculations in which both possibilities are entertained in order to consider their impact on 
thermodynamics and their astrophysical implications. 



C. Finite temperature case 

We now investigate results at finite temperature and values of the lepton content characteristic of those 
likely to be encountered in the evolution of a PNS. We choose three representative sets of PNS conditions 
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which correspond to: the initial conditions within a PNS (entropy/baryon s = 1, trapped neutrinos with 
a lepton fraction = 0.35), a time after several seconds when the interior is maximally heated (s = 2, no 
trapped neutrinos so = 0), and a very late time when the PNS has cooled (s = 0, F,^ = - identical to 
the zero temperature case discussed above). For a detailed explanation of the evolution of a cooling PNS 
see Pons et al. [Q. 

The contribution of the nucleons to the entropy per baryon sj^ = Sj\f/nB, with = + Up denoting 
the total nucleon density, in degenerate situations [T/Ep^ ^ 1) can be written 



SAT = vr T ^ ,3 , (49) 

where M* and kp^i are the effective mass and the Fermi momentum of species i, respectively. For the 
temperatures of interest here, and particularly with increasing density, the above relation provides an 
accurate representation of the exact results for entropies per baryon even up to = s„ + Sp ~ 2. The 
behavior with density of both the Fermi momenta and the effective mass controls the temperatures for a 
fixed sj^. 

For kaons it is straightforward to show that the contribution to the entropy from mesons can be 
ignored since it is exponentially suppressed in comparison to the K~ contribution. For the latter, keeping 
the leading temperature dependence of the simplest approximation scheme for bosons given in Ref. []36| , 
the kaon entropy per baryon is 
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and y is determined from ip by solving the equation 

-0 = 1 — y + lny. (51) 

Below the kaon condensation threshold as the temperature becomes very small t/^^— oosoy— s>0. Above 
the kaon condensation threshold the last of Eqs. ( |30D implies that t/^ = in which case y = 1. This 
simple approximation provides quite an accurate account of the kaon entropy per baryon which is fairly 
small for the scenarios examined here since it involves just the thermal contribution and the condensate 
plays no role. The total entropy per baryon Stot = sj^ + sk + [Se + + Su)/nB also includes the lepton 
contributions; Stot is dominated, however, by sj\f. 

In Figs. and |12|, the relative concentrations of various particles are displayed versus baryon number 



density for our three PNS conditions in the cases GM+GS and ZM+TW, respectively. The cases shown 
allow the Gibbs equations to be solved, and the boundaries of the mixed phase regions are indicated by 
vertical lines. The effect of finite temperature is to allow the existence of /i~ and K~ particles at all 
densities, although kaons become relatively abundant only within the mixed phase region. In the third 
set of diagrams, trapped neutrinos are present at all densities and the appearance and abundances of 
the negatively charged particles /i~ and K~ are suppressed. Furthermore, the critical density for kaon 
condensation is shifted to higher density. 



In Fig. the pressure is displayed as a function of baryon number density for these two Lagrangians 
and the three PNS conditions. Two choices for the kaon optical potential are shown to highlight differences 
between cases in which kaons condense in second or first order phase transitions. The reduction of the 
pressure when kaons condense is obvious. For conditions in which the phase transition is first order, the 
result of applying the Gibbs conditions and the result of assuming pure phases (thin line) are both shown. 
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The application of the Gibbs conditions leads to further softening of the pressure over a wider density 
range. In the case of model ZM+TW, a first order phase transition occurs only for very low temperatures 
and low neutrino concentrations. 



In Fig. |14| we show the matter temperature as a function of the baryon density for these two La- 



grangians for the two PNS conditions with s > (the kaon optical potentials are as in the previous figure). 
The appearance of the kaon condensate generally leads to a reduction in specific heat which is indicated 
by the abrupt temperature increase which persists to high densities. In the case of first order transitions, 
applying the Gibbs conditions leads to a further enhancement of the temperature in the mixed phase 
regime. This behavior is in marked contrast to the case in which additional fermionic degrees of freedom, 
such as hyperons or quarks, are excited [^] causing the temperature to drop and the specific heat of the 
matter to be increased. The latter follows from Eq. (^) where, in the absence of any variation of M*, a 
system with more components at a given baryon density has a smaller temperature than a system with 
fewer components (recall that J^i^i = 1). In the present case the dropping of the effective mass is the 
dominant effect and this leads to larger temperatures. 

Figure ^ shows the phase diagram of kaon condensed matter, for the case GM+GS with = —120 
MeV. The left panel displays results for zero temperature in the density-lepton concentration plane. The 
dashed lines show the minimum lepton concentration allowed at zero temperature (with Yi^ = 0) for each 
density. Note that the minimum lepton concentration increases with density until the phase transition 
begins; above this density, the minimum lepton concentration decreases with increasing density. Also note 
that the phase transition to a kaon-condensed phase is pushed to higher densities when trapped neutrinos 
are present. This implies that in the initial PNS core material, in which Yl ^ 0.35 — 0.4 and the central 
density is less than 3.5 times the nuclear saturation density, a kaon condensate phase likely does not 
exist. However, as neutrinos leak from the star the transition density decreases and a kaon condensate 
eventually forms. The right panel displays results in the density versus temperature plane, assuming no 
trapped neutrinos {Y^, = 0). 

The phase diagram for kaon condensed matter for the case ZM+TW with U^j^ = —140 MeV is shown 
in Fig. |T^; the results are qualitatively similar to the GM+GS case in which U]!^ = —120 MeV in Fig. 
|15|. This is understandable from the perspective that the actual optical potential for these two models 
are nearly the same. The boundary between phases I and the mixed-phase region are nearly the same 
for the two cases. The major difference is the much smaller extent of the mixed-phase region for the case 
ZM+TW. 

Note that for both cases the density at which the phase transition begins is relatively independent 
of temperature, so that the heating which initially occurs in the PNS has little effect on the eventual 
appearance of a kaon condensate. Also note that the density range of the mixed phase decreases with 
increasing temperature, and the mixed phase persists to high temperatures. It appears that the mixed 
phase exists up to temperatures exceeding 60 MeV, for the case GM+GS and = —120 MeV, or 30 
MeV for the case ZM+TW with f/)™ = —140 MeV. It becomes increasingly difficult to determine the 
properties of a mixed phase near the temperature at which it disappears. 



In Fig. O the gravitational mass is plotted as a function of central baryon number density for these 



models. Results are shown for our three PNS conditions which correspond to snapshots of the PNS 
evolution. The initial configuration (dotted curves) has the largest maximum mass. The progression 
to the dashed and solid curves indicates the evolution with time and we see that the maximum masses 
decrease. The effect of temperature upon the structure of the PNS is significant. Thermal kaons play 
a significant role here, since the net negative charge they contribute to the system partially inhibits the 
appearance of the condensate which allows hot neutrino-free stars to reach higher masses than cold stars. 
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The net decrease in maximum mass during the evolution for either case is seen to be of order 0.2-0.3 M©. 
Thus there is an appreciable range of masses for the PNS which will result in metastability with the star 
ultimately collapsing to a black hole. The central density of the maximum mass, zero temperature star 
is smaller for the GM+GS case than for the ZM+TW case. This is in spite of the apparently "softer" 
GM+GS equation of state in which the kaon condensed mixed-phase region extends over a wider density 
range. Ultimately, the smaller maximum mass of the GM+GS EOS leads to a smaller central density at 
the maximum mass. 

IV. SUMMARY AND OUTLOOK 

In this work, we have studied the equation of state of matter, incorporating the possible presence of a 
kaon condensate, and including the effects of trapped neutrinos and finite temperatures. The calculation 
of the neutrino spectra of different flavors emitted from a proto- neutron star as it evolves from a hot, 
lepton-rich state to a cold, neutrino-poor state requires the knowledge of the equation of state of matter at 
temperatures up to about 50-60 MeV and lepton fraction up to about 0.4. Since the nucleon-nucleon and 
kaon-nucleon interactions at high density are relatively poorly understood, we explored several possible 
field-theoretical models in both sectors. These models are distinguished by the form of the assumed scalar 
(and in some cases vector) interactions which chiefly determine the density dependences of the miclcon 
and kaon effective masses. These models produce significantly different high density behavior of the EOS, 
even though the kaon-meson couplings in these models are calibrated to give the same the kaon-nucleus 
optical potential in nuclear matter. 

The principal findings of our studies at zero temperature were: 

1. The order of the phase transition between pure nucleonic matter and a phase containing a kaon 
condensate depends sensitively on the choice of the kaon-nucleon interaction. 

2. In one case we studied (KPE), although a first-order phase transition resulted, it was not possible to 
satisfy Gibbs' rules for phase equilibrium which would have produced a mixed phase. We performed 
a detailed analysis of this situation and found that scalar, and to a lesser extent the isovector, 
interactions that vary rapidly with density were chiefly responsible for this failure. This was conflrmed 
by developing a new kaon-nucleon interaction (TW) with more moderate variations in the scalar 
density and the kaon effective mass in which the Gibbs' criteria in a first order phase transition 
would be satisfied. The extent of the mixed phase region was thereby reduced. The significance of 
the new kaon-nucleon interaction (TW) we developed is that it avoids the anomalous behavior for 
the kaon effective mass that occurs in previous models (KPE, GS) at very high density. Near the 
low-density boundary of a mixed phase region, the kaon condensed phase appears with large density, 
too large for the KPE interaction to produce physically acceptable effective masses. We also made 
detailed comparisons with earlier work which used the GS form for the scalar interactions. 

3. In all models considered (KPE, GS and TW), a first-order phase transition occurs only for large 
values of the kaon-nucleus optical potential; moderate values generally produce a second order phase 
transition. 

In the meson exchange models studied here, only linear kaon self-interactions were considered. In the 
case of a first order phase transition, the condensate amplitude was found to be rather large at the low- 
density boundary of the mixed phase. We therefore explored the effect of non-linear kaon self-interactions 
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guided by the chiral model in Appendix B. We found that introducing higher order interactions, using 
the lowest order chiral Lagrangian, results in a second order, rather than a first order, phase transition. 
Whether this behavior persists when more general higher order operators in the chiral expansion are 
considered remains an open question. 

At finite temperatures, we find the effects of condensation, in general, are less pronounced than at 
zero temperature. For moderate values of the optical potential, when the phase transition is first order at 
zero temperature, kaon condensation eventually becomes a second order phase transition at high enough 
temperatures, whether or not neutrino trapping is considered. The temperature at which this occurs is 
in the range of 30-60 MeV, depending upon interactions. For the cases at finite temperatures in which 
the transition is first order, its thermodynamics (such as the pressure- density relation) becomes effectively 
similar to that of a second order phase transition. This is because of the existence of thermal kaons and 
because of nucleonic thermal effects. The condensate is suppressed, and moved to higher densities, both 
by the existence of trapped neutrinos and by finite temperatures. Compared to earlier works, the new 
aspects of our work are: 

1. The delineation of the phase boundaries in the baryon density versus lepton number and baryon 
density versus temperature planes. This is helpful to anticipating the possible outcome in a full PNS 
simulation. In particular, the critical temperatures above which the mixed phase disappears are 
above 30 and 60 MeV, depending upon the interaction. This has implications for the temperature 
dependence of the surface energies, and for the melting temperatures, of the droplets in the mixed 
phase. 

2. The finding that thermal effects on the maximum gravitational mass of neutron stars are comparable 
to the effects induced by the trapped neutrino content. This is in stark contrast to previously studied 
cases in which nucleons-only matter, or matter containing hyperons, were considered. Furthermore, 
compared to equations of state previously studied that allow metastable protoneutron stars, those 
containing hyperons or quark matter, the maximum mass does not significantly decrease during the 
deleptonization of the protoneutron star because of these thermal effects. Only after the temperature 
in the protoneutron star significantly decreases does the maximum mass appreciably fall. This implies 
that the possible collapse of a metastable protoneutron star to a black hole occurs during the late 
stages of cooling, after several tens of seconds, rather than during the late stages of deleptonization, 
which is somewhat earlier. 
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APPENDIX A: MESON EXCHANGE VERSUS CHIRAL MODELS 

In this Appendix, we examine the conditions under which there exists a close correspondence between 
a meson exchange model and the chiral SU{3)l x SU{3)ji approach of Kaplan and Nelson p. Such a 
correspondence is most easily established for the zero temperature case by setting the scalar self-coupling 
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terms, i.e., U{(j) = 0. Specializing to the case where the only baryons are nucleons and using the Walecka 
Lagrangian for the nucleons, it was shown in Ref. ^ that the chiral thermodynamic potential per unit 
volume can be written 



rij^ + 12/2 12/2 1 2j/2 I o f 
y = l^a^^ - 2"^-^0 - 2"^P^ +^2^]j2 



+2mlf sin'le - Wf sin^ 9 , (52) 



where the primes on the meson fields distinguish them from those used previously and B is the Heaviside 
step function. The nucleon effective masses are 

M* = M- g^a' + (2a2 + 4a3)m, sin^i^ 

M* = M- gy + (2ai + 2a2 + 4a3)m, sin^i^ . (53) 

We employ the values suggested by Politzer and Weise [|10|, namely ainis = —67 MeV (rris is the strange 
quark mass) and a2ms = +134 MeV. a^rris is usually taken to lie in the range —134 to —310 MeV. If we 
ignore the fairly small effect of aim^ here and in the kaon-nucleon sigma term, = — i(ai+2a2+4a3)ms, 
we can write 

M* c:^ M; c:^ M* = M - g,a' - 2S^^ sin^i^ ^ m - g„a 

E^p ~ = VWTm^ . (54) 
As well as redefining the scalar field, we can redefine the chiral vector fields entering the chemical potentials: 

/ip = t'p + gw^'o + \gph'Q - 2/x sin^|6' = z/p + g^ujQ + \gphQ . (55) 



Substituting in Eq. ( ^21) we find 

Vt}j- -\- VLk 1 2 2 1 2 2 1 27,2 , o /" /TP* \r\r n^* 




^'C-^^'^]}^ (56) 



If we expand this in powers of 6 and retain only the lowest order 6"^ term, the last term in Eq. (|5^) does 
not contribute and our thermodynamic potential is exactly of the form given by Eqs. (^) and (p6| ) for the 
meson exchange model provided that the axpE expression is used. In order for the correspondence to be 
exact, the parameters for the a and u meson need to obey 

9a9aK _ _ gu)gu)K _ _3_ ^^^^ 

These are precisely the conditions found in Ref. |^ for the optical potentials of the chiral and meson 
exchange models to be the same in nuclear matter. The relation involving the uj meson couplings is quite 
well obeyed with our parameters. In addition, for the p meson. 



20 



^ = 77i • (58) 

This indicates that {gpQpK) I {Qi^gLuK) — f; ^ condition which is not well obeyed by the parameters used 
here or in other works. Given that the chiral and meson exchange thermodynamic potentials can be 
put into precise correspondence to lowest order in 6'^, it follows that the equations of motion and the 
thermodynamics will be identical to this order. 

If scalar self-coupling terms are included, U{a) ^ 0, then the transition from the chiral to the meson 
exchange approach will couple higher powers of the a field to the kaon condensate (in the braces in Eq. 
(^)). It will also introduce higher order terms. These additional contributions may not be negligible so 
the correspondence between the two approaches becomes less precise. 



APPENDIX B: HIGHER ORDER KAON SELF-INTERACTIONS 

Our findings in Appendix A naturally raise the question of whether it is sufficient to work at order 
involving only linear kaon self-interactions, in the meson exchange models. It clearly will be sufficient at 
the low-density onset of a second order phase transition where 9 is small. On the other hand, for a first 
order phase transition, the value of 9 is large at the low-density onset of the mixed phase, particularly for 
the GS model. It is therefore interesting to explore the effect of non-linear kaon self-interactions guided 
by the chiral model. 

The order of the phase transition (in the mean field approximation) is determined by the behavior of 
the thermodynamic potential, Vl{9), at fixed chemical potentials. A first order transition, with a mixed 
phase, is possible only if there exists some value of yU„ for which Vt{9) exhibits two degenerate minima. At 
the critical density corresponding to the low-density onset of the mixed phase, the 6' = phase should be 
a local minimum which is degenerate with a minimum at some finite 9 = 9^. In the vicinity of the critical 
density, the 9 = phase is nearly charge neutral (with an infinitesimal excess of positive charge and a 
volume fraction close to 1 which balances the negative charge in the kaon phase which has an infinitesimal 
volume fraction). This requirement enables us to determine the electron chemical potential at the critical 
density by charge neutrality. 

We focus on the GM+GS model for which the thermodynamic potential of nucleons and kaons was 
given in Eqs. (|^) and (|26|) ; the contribution due to leptons is ignored since it does not contain any 9 
dependence at fixed fi. At zero temperature with a kaon optical potential = — 120 MeV, this model 
predicts a first order phase transition in the vicinity of fin = fJ'c — 1160 MeV, as can be deduced from 
Fig. |T^ where Q{9) — fl{9 = 0) is plotted as a function of 9. The thermodynamic potential for the model 
GM-I-GS is shown as the solid curve labelled It clearly shows two minima, one at 6' = and the other 
at ^ ~ 2. The latter corresponds to a kaon number density ~ 1 fm~^ which is larger than the baryon 
density. For such a dense condensate one would suspect that non-linear kaon self-interactions might be 
important. The order 9'^ corrections to the thermodynamic potential are easily found from Eq. (|56|) : 



Afi4 = -^«-V). (59) 



The result of adding this correction to is shown as the dashed curve in Fig. [Tg. It greatly alters the 
behavior of Q{9) for 9> 1. The exsistence of a second minimum suggests that a first order phase transition 
is still possible, but at larger However, we find that this is not the case and a second-order transition 
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occurs at /i„ = 1213 MeV. It is possible to incorporate all powers of 6 arising from self-interactions in the 
chiral model. In this case the correction to the grand potential is 

AQn = 2ml f sin^l^ - i/iV' sin^ 6 - - /i^) . (60) 



The result of including this correction is shown as the dot-dashed curve in Fig. 0. In this case no first 
order phase transition is possible in the vicinity of /i„ = Hc- Instead a second order phase transition 
occurs once again at = 1213 MeV; this is because kaon self interactions play no role when 9 is small. 
Despite our findings here, it is not clear if kaon self-interactions will generically disfavor a first order 
transition. This is because we have ignored higher order operators in the chiral expansion which will 
become important with increasing 9. The indication from phenomenological chiral perturbation theory 
p7| is that such effects can be significant when 6' ~ 2. The robust finding here is that the higher order 
kaon self-interactions predicted by the lowest order chiral Lagrangian lead to a second order, rather than 
a first order, phase transition. 
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TABLE I. Meson-nucleon coupling constants fitted to a binding energy /particle of 16.3 MeV at an equilibrium 
density of Hq = 0.153 fm~^ in nuclear matter with a compression modulus K = 240 MeV and effective mass 
M* = 0.78M. The symmetry energy coefficient agym is 32.5 MeV. 



Model 


9a 
rUa 

(fm) 


Qui 

(fm) 


9p 

nip 

(fm) 


b 


c 


GM 


3.1507 


2.1954 


2.1888 


0.008659 


-0.002421 


ZM 


3.1228 


2.1954 


2.1888 


-0.006418 


0.002968 



TABLE II. Kaon optical potentials for the models GS and TW. The values in the second column refers to 
the linear approximation of the exact results for the models GS and TW shown in the third and fourth columns. 
All results in this paper have been labelled according to the linear approximation (in order to make a comparison 
with the results of GS), which is exact for the model KPE. 



gaKO- 


rrlin ttKPE 
— Uj^ — —Uj( 




ttTW 


63 


80 


81 


80 


103 


100 


100 


95 


143 


120 


117 


109 


183 


140 


136 


122 


223 


160 


150 


134 
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FIGURE CAPTIONS 



FIG. |I|: Pressure versus baryon number density for the six choices of the nucleon and kaon Lagrangians 
considered in this paper. The temperature T = and there are no trapped neutrinos (Yi, = 0). Selected 
values for the kaon optical potential are indicated. The solid lines show the pressure in the pure 
phases I (nucleons only) and II (the high-density nucleon- kaon condensed phase). The dashed lines show 
the pressures obtained by imposing Gibbs' criteria for phase equilibrium in a mixed-phase region for the 
case of first order transitions. For the KPE choice of the kaon Lagrangian, Gibbs' criteria could not be 
satisfied despite the occurence of first order phase transitions in some cases. 

FIG. ^ The density dependences of the scalar, vector, and iso-vector fields for different choices of the 
nucleon and kaon Lagrangians (T = 0,Y^ = 0). The solid curves show the chemical potential fin — fJ'p = f^- 
In this figure, results are shown only for the pure phases I and II; the mixed phase produced by satisfying 
Gibbs' criteria is ignored. 

FIG. ||: As for Fig. ||, but for the density dependences of the kaon and nucleon effective masses. The 
solid curves show the condensate amplitude. 

FIG. The electron chemical potential fi versus the neutron chemical potential /i„ in pure phases 
I and II, and in the mixed phase. The pure phase I (solid curve) consists of nucleons and leptons. The 
pure phase II (dashed curve) is comprised of a kaon condensate coexisting with nucleons and leptons. The 
mixed phase (dots) is constructed by satisfying Gibbs' rules for phase equilibrium. 

FIG. If Individual charge densities of pure phases I and II and the volume fraction x of phase I in the 
mixed phase as a function of baryon density. Results are for the GM+GS model with = —120 MeV 
and for ZM+TW model with U^^"" = -140 MeV. 

FIG. 1^: The density dependences of the scalar, vector, and iso-vector fields for two choices of the 
nucleon and kaon Lagrangians (T = 0,Yi, = 0). Phase I is the pure nucleon phase and phase II is the 
high-density nucleon-kaon condensed phase. The vertical lines demark the mixed phase region. 

FIG. ^ The density dependences of the chemical potential fin — fJ^p = fi, the kaon (K) and nucleon 
(N) effective masses, and the condensate amplitude (T = 0, Fj, = 0). Notation is as in Fig. ^j. 

FIG. |: The electron chemical potential fi in phase II matter versus charge density for different models 
at a fixed neutron chemical potential of /i„ = 1250 MeV. In all cases, the optical potential = —120 
MeV. 

FIG. ^ Left panel: The gravitational mass as a function of the central baryon number density for 
the cases GM+GS and ZM+TW (T = 0, = 0). Curves are labelled by the values of ?7f^ and the EOS 
includes a mixed phase region. Right panel: The gravitational mass as a function of the stellar radius. 

FIG |10|: The density dependences of the kaon energy u in matter with different isospin content. The 
top panel compares results of GM+KPE, GM+GS and ZM+TW models for beta stable neutron star 
matter for = —80 and -120 MeV, respectively. The bottom panel shows results for the ZM+TW 
model with = -80 MeV in pure neutron matter, beta stable neutron star matter and nuclear matter. 

FIG 0: The relative concentrations of hadrons and leptons as functions of baryon number density 
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for three representative snapshots during the evolution of a PNS. The results shown are for the model 
GM+GS. To the left of the vertical line there is no kaon condensate, to the right a mixed phase is present. 



FIG pi Same as Fig. O, but for the model ZM+TW. 



FIG 13: The pressure versus baryon number density for three representative snapshots during the 



evolution of a PNS. The cases shown in the upper panels produce only second order phase transitions. For 
the cases in the lower panels the transitions are first order, except for ZM+TW with s > 0. In the lower 
panels, heavy curves include a mixed phase region and light curves ignore a mixed phase region. 

FIG The temperature as a function of baryon density for two snapshots during the PNS evolution. 



Other notation is as in Fig. [12. 



FIG |T5|: The phase diagram of kaon condensed matter for the case GM+GS and = —120 MeV. 
The left panel shows results at zero temperature in the density versus lepton concentration plane. The 
dashed curve shows the minimum lepton concentration for each density, which occurs for trapped neutrino 
concentration Y^, = 0. The right panel shows results in the density versus temperature plane for neutrino 
free matter (Fj, = 0). 



FIG ra Same as Fig. 15, but for the model ZM+TW and 



lin 
K 



-140 MeV. 



FIG The gravitational mass versus central baryon number density in the GM+GS and ZM+TW 
models for three representative snapshots during the PNS evolution. 



FIG 0. The thermodynamic potential as a function of the condensate order parameter 9. Results are 
shown for the GM+GS model near the critical density (/i^ ~ 1160 MeV and ~ 243 MeV) with optical 
potential f/f^ = -120 MeV. 
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FIG. 3. 
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